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REPORT  ON  A  SEMINAR  ON  THE  HYDRODYNAMIC  THEORY  ASSOCIATED 
WITH  SHIP  MOTION  IN  WAVES 
by 

Paul  Kaplan  and  Jack  Kotik 


1.  INTRODUCTION 

The  theory  of  the  motions  of  a  ship  in  waves  was  developed 
along  two  major  diverse  paths.  Various  practical  investigators  and 
some  designers  have  used  results  based  on  many  questionable  assumptions, 
while  more  academic  investigators  have  studies  isolated  parts  of  an 
overall  motion  analysis  (e.g.  damping  and  added  mass  coefficients)  or 
have  treated  idealized  problems  having  no  direct  application  to 
practice.  As  an  example  of  practical  approaches  the  work  of  Korvin- 
Kroukovsky  and  Jacobs^  ^  develops  a  set  of  coupled  linear  differential 
equations  for  heave  and  pitch  of  a  ship  in  regular  head  seas,  and  the 
coefficients  used  there  are  the  values  obtained  by  "strip  theory",  i.e. 
by  use  of  integrals  of  two-dimensional  quantities.  A  number  of  different 
methods  have  been  used  to  calculate  the  two-dimensional  coefficients 
which  appear  in  such  sets  of  differential  equations,  each  having 
different  degrees  of  validity  and  simplicity.  Similarly  there  exists 
a  number  of  different  expressions  for  hydrodynamic  forces  on  oscillating 
bodies  advancing  on  or  below  the  free  surface,  which  have  been  obtained 
within  what  is  considered  to  be  a  three-dimensional  analysis  (subject 
to  certain  simplifications  such  as  thin- ship,  slender-body,  symmetry, 
etc.) 

It  seemed  appropriate  to  bring  together  a  group  of  investiga¬ 
tors  who  were  active  in  the  various  theoretical  aspects  of  the  problem, 
and  who  represented  some  of  the  different  approaches  to  the  problems 
presently  being  investigated,  in  order  to  conduct  an  informal  seminar 
on  this  subject.  The  participants  in  the  seminar  and  their  affiliations 
are  listed  below: 


TECHNICAL  RESEARCH  GROUP 


2. 


Dr.  Klaus  Eggers,  Institut  fur  Schiffbau,  Hamburg  University 

Dr.  Paul  Kaplan,  TRG,  Inc.* 

Dr.  Samuel  Karp,  Hew  York  University;  consultant  to  TRG,  Inc. 

Dr.  Jack  Kotik,  TRG,  Inc. 

Dr.  Richard  C.  MacCamy,  Carnegie  Inst,  of  Tech. 

Dr.  J.  Nicholas  Newman,  David  Taylor  Model  Basin 

Dr.  Reiner  Tinman,  Technological  University,  Delft 

These  scientists  were,  in  the  main,  primarily  concerned  with  problems 
related  to  the  theoretical  calculation  of  various  hydrodynamic 
coefficients  used  in  the  equations  of  motion,  rather  than  with  the 
complete  formulation  of  equations  for  calculation  of  actual  motions. 
While  certain  of  the  participants  had  experience  in  the  latter  area, 
the  majority  interest  led  to  concentration  on  the  more  theoretical 
hydrodynamic  aspects  of  the  ship  motion  problem.  However,  some  atten¬ 
tion  was  paid  to  the  formulation  of  different  types  of  equations  of 
motion,  and  other  aspects  closer  to  practical  utility.  Nevertheless, 
it  still  remains  to  integrate  the  results  of  this  seminar,  and  work 
that  was  produced  or  will  be  produced  as  a  result  of  its  recommendations, 
into  more  effective  methods  for  determining  the  motion  of  a  ship  in 
waves. 

The  general  topics  covered  were  roughly  as  indicated  below: 

a)  Damping  and  inertia  coefficients  for  the  forced  periodic 
motion  of  a  body  at  zero  speed. 

b)  Influence  of  forward  speed  on  the  above. 

c)  Cross-coupling  terms  between  different  degrees  of 
freedom. 

d)  Evaluation  of  forcing  functions  due  to  waves  incident 
on  a  moving  ship. 

e)  Equations  of  motion,  including  convolution-type 
integro-differentlal  equations. 


Currently  at  Oceanic 8,  Inc. 
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The  main  emphasis  was  on  heave  and  pitch  motions,  but  some  considera¬ 
tion  was  given  to  rolling  and  the  general  problem  of  lateral  motions. 

No  attempt  was  made,  either  in  the  seminar  or  in  this  report,  to  arrive 
at  an  encyclopedic  summary  of  past  work  and  the  present  state  of 
knowledge  .  The  participants  felt  free  to  devote  their  attention  to 
problems  which  most  interested  them,  and  as  a  result  some  important 
questions  were  considered  briefly  or  not  at  all.  In  the  case  of  those 
problems  which  were  discussed  at  length  considerable  attention  was 
paid  to  the  practical  utility  of  the  existing  and  anticipated  results. 

This  seminar  was  conducted  un  er  the  sponsorship  of  the 
Office  of  Naval  Research  under  Contract  Nonr-3175(00) . 


Vossers^  has  published  an  extensive  summary  of  the  subject  which 
can  be  recommended  as  a  starting  point. 
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2.  SUMMARY  OF  DISCUSSIONS 

2.1  Present  Analytical  Methods  for  Predicting  the  Motions  of 

Ships  In  Waves 

A  simple  method  of  analysis  for  the  problem  of  a  ship  in 
regular  waves  is  to  assume  that  it  is  a  linear  uncoupled  dynamic  system 
in  its  various  degrees  of  freedom.  Since  such  equations  are  easily 
solved  the  main  problem  is  to  determine  the  appropriate  coefficients, 
including  the  oscillatory  forcing  functions,  from  available  hydrodynamic 
theory.  An  example  of  this  approach  is  the  paper  by  Weinblum  and 
St.  Denis ^  where  the  hydrodynamic  techniques  used  were  the  Froude- 
Kryloff  hypothesis  (the  assumption  that  the  forces  due  to  waves  are 
the  same  as  those  exerted  on  the  water  displaced  by  the  ship,  i.e. 
neglecting  interactions  between  fluid  and  ship),  two-dimensional  damping 
and  added  masses  for  ship  sections  which  were  integrated  over  the  ship 
length  in  accordance  with  strip  theory,  and  the  use  of  infinite- fluid 
finite- length  inertia  factors  for  three-dimensional  corrections.  No 
comparisons  with  experiment  were  made  for  any  practical  ship  form  in  [2]. 

A  more  complicated  set  of  linear  equations,  including  coupling 
terms,  showing  fair  agreement  with  experiment  for  five  ship  forms,  was 
given  by  Korvin-Kroukovsky  and  Jacobs^  .  The  technique  used  there  is 
appealing,  but  the  evaluation  of  the  various  hydrodynamic  coefficients 
used  in  the  equations  is  still  the  subject  of  much  discussion  and 
research.  In  addition*  certain  questions  arise  as  to  the  adequacy  of 
the  equations  themselves  when  considering  non- monochromatic  waves  as 
well  as  problems  of  nonlinearity.  While  these  equations  have  been 
successful  in  certain  cases,  and  of  course  only  for  problems  of  heave 
and  pitch  in  monochromatic  head  seas,  the  need  for  appropriate  tools 
with  wider  applicability  with  regard  to  ship  form  and  number  of  degrees 
of  freedom  is  certainly  evident.  However,  it  is  desirable  to  extend 
their  equations  to  six  degrees  of  freedom  and  monochromatic  waves  of 
arbitrary  (single)  heading,  and  to  obtain  more  accurate  coefficients. 

This  accounts  for  the  great  Interest  at  the  seminar  in  the  evaluation 
of  damping  and  inertia  coefficients,  coupling  terms,  forcing  functions, 
etc. 
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2.2  Present  Methods  of  Determining  the  Damping  and  Inertia 
Coefficients  In  the  Equations  of  Motion  of  a  Ship  In 
Regular  (monodirectional. monochromatic)  Waves. 

2.2.1  Two-Dimensional  Techniques 

Considerable  effort  has  been  devoted  to  the  two- 
dimensional  problem  of  the  forces  acting  on  a  section  oscillating  at 
the  free  surface.  The  results  are  expressions  for  the  damping  coefficient 
(or  wave  amplitude  coefficient)  and  the  added  mass  coefficient  of  the 
section,  as  functions  of  a  dimensionless  frequency  parameter.  The 
various  techniques  used  to  determine  these  quantities  are  described  in 
the  following  brief  outline. 

A  simple  method  used  for  determining  the  damping  due 
to  heaving  motion  of  a  section  on  the  free  surface  is  due  to  Havelock 
and  Holstein,  ^  and  makes  use  of  a  distribution  of  sources  along 
the  body  surface,  with  the  strength  proportional  to  the  normal  velocity 
along  the  body  contour.  The  radiated  wave  amplitude  is  found  and  this 
leads  to  a  value  of  the  damping  coefficient.  A  more  useful  and  also 
more  valid  method  of  obtaining  both  sectional  damping  and  added  mass 
is  that  of  Grim^,  which  is  applied  to  a  mathematical  family  of 
sections  known  as  the  Lewis  forms.  Grim  satisfies  the  boundary  condi¬ 
tion  on  the  body  at  a  finite  number  of  points.  He  has  published  the 
results  for  a  variety  of  sections.  Information  for  other  sections 
can  be  obtained  by  interpolation  from  Grim's  charts. 

One  of  the  first  rigorous  treatments  of  a  linearized 
forced  motion  problem  was  Ursell's  treatment  ^  of  the  heaving  semi- 
immersed  circular  cylinder.  The  problem  was  transformed  to  that  of 
solving  an  infinite  linear  system.  Convergence  was  proved,  and  the 
wave  amplitude  was  computed  numerically  by  solving  a  truncated  set 
of  linear  equations.  More  recently  Porter^7 ^  has  extended  Ursell's 
method.  He  transformed  the  problem  of  a  heaving  section  to  that  of 
solving  an  infinite  linear  system  whose  coefficients  Involve  the 
coefficients  in  the  series  expansion  of  the  conformal  mapping  of  the 
exterior  of  the  section  on  the  exterior  of  a  circular  section.  He 
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programmed  the  IBM  704  computer  to  solve  truncated  sets  of  linear 
equations,  and  computes  wave  amplitude,  virtual  mass  and  the  pressure 
on  the  section  at  discrete  points.  The  number  of  equations  treated 
is  an  input  variable  and  by  increasing  the  number  until  the  results 
converge  one  can  be  assured  of  a  quite  accurate  solution.  In  the  case 
of  a  circle  the  agreement  with  Ursell  is  excellent;  which  is  to  be 
construed  as  a  verification  of  Ursell' s  results.  Results  for  a 
number  of  other  sections  are  presented  in  [7],  and  the  computer  program 
can  be  used  for  additional  cases  as  required.  In  conjunction  with  a 
program  for  determining  conformal  mappings  numerically,  the  Porter 
program  could  be  used  to  treat  arbitrary  sections. 

Kaplan  and  Jacobs^  computed  added- mass  and  damping 
coefficients  for  sections  by  a  perturbation  method  assuming  small  beam. 
While  the  work  does  show  good  agreement  with  known  results  such  as  those 

for  a  circular  or  elliptical  section,  it  is  known  that  there  is 
an  error  of  the  order  of  the  vertical  added  mass  coefficient  in  the 
asymptotic  value  at  large  frequencies. 

Tasai1  4  has  computed  virtual  mass  and  damping  for 
a  number  of  simple  forms  by  a  method  which  appears  to  resemble  Porter's 
except  that  the  calculations  were  not  mechanized.  Hfs  results  are 
generally  very  close  to  Porter's,  although  there  are  some  deviations. 

MacCamy^^  has  developed  a  perturbation  theory  for 
sections  of  small  draft.  The  integral  equation  to  which  he  reduced 
the  problem  has  been  treated  numerically  for  a  heaving  strip.  A 
third  solution  was  obtained  by  Porter  (see  Figure  1.)  In  addition, 

MacCamy  has  investigated  related  problems  such  as  the  behaviour  of  the 
potential  at  the  confluence  of  two  different  boundary  conditions,  and 
the  behaviour  of  the  potential  at  low  frequencies.  To  appreciate  the 
second  problem  note  that  the  problem  at  zero  frequency  is  a  Neumann 
problem  whose  solution  is  not  unique.  Yet  for  every  non-zero  frequency 
the  potential  is  uniquely  determined.  To  which  of  the  many  solutions 
of  the  zero- frequency  problem  does  the  potential  tend  a a  the  frequency 
tends  toward  zero  ?  MacCamy  has  sketched  a  proof  that 
$  -  0Q  +  (of41ogK)  /  Vnds  +  O(KlogK)  as  k  -*-  0,  where 
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FIGURE  I 

RATIO  OF  WAVEHEIGHT  TO 
HEAVE  AMPLITUDE  FOR  A 
STRIP.  PORTER'S  RESULT 
IS  EXACT  EXCEPT  FOR  NU¬ 
MERICAL  TRUNCATION 
ERRORS. 
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K  -  c2/g  -  2 ir/X 

0Q  -  any  zero- frequency  potential 
a  ■  a  constant  that  depends  on  the  choice  of  0Q 
Vn  ■  prescribed  normal  velocity  on  the  section. 

For  a  symmetrical  motion  (heave)  of  a  symmetrical  section  f  V^ds  y  0 
in  general  so  that  an  "infinite  constant"  appears.  This  accounts  for 
the  logarithmic  infinity  in  the  virtual  mass  of  a  heaving  section  at 
zero  frequency,  but  it  does  not  affect  the  leading  term  in  the  damping 
coefficient  at  low  frequency. 

The  asymptotic  behaviour  of  the  wave  amplitude  and 
added- mass  coefficients  at  low  frequency  has  been  investigated  by  a 
number  of  authors  with  varying  degrees  of  rigour.  It  has  been  known 
for  some  time  that  the  wave- amplitude  coefficient  for  the  heaving  motion 
of  an  arbitrary  section  is 

2Ka, 

where  2a  is  the  width  at  the  waterline.  Ursell^^  has  shown  hot;  to 
find  the  leading  term  in  the  wave  amplitude  for  the  heaving  and  rolling 
motion  of  sections  for  which  the  mapping  function  is  known.  He  also 
found  the  leading  term  of  the  added-mass  coefficient  for  heave  at  low 
frequency. 

The  high-frequency  asymptotics  has  been  investigated 
by  Ursell  by  transforming  the  usual  integral  equation  of  the  problem 
into  another  integral  equation  whose  kernel  is  small  at  high  frequency. 
The  method  was  applied  to  a  heaving  circular  cylinder  and  a  heaving 
sphere,  with  the  following  results: 

Cylinder: 

Wave  amplitude  coefficient:  ^ 

Added- Mas 8  coefficient:  1  -  3^5 

TRG  has  proposed  a  simple  formula  for  the  leading  term  of  the  wave- 
amplitude  coefficient  at  high  frequency  in  terms  of  the  local  behaviour 
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of  the  infinite-frequency  potential  at  the  intersection  of  the  section 
and  the  free  surface.  The  method  has  been  applied  to  the  following 
cases: 

Elliptic  Cylinder:  Roll  and  Heave 
Horizontal  Strip:  Roll  and  Heave 
Vertical  Strip:  Roll 
Lenticular  Cylinder:  Heave 

The  results  agree  with  Ursell  in  all  cases  for  which  comparison  is 
possible. 

A  substantial  forward  step  is  the  introduction  of  the 

Hi-Fi  method  for  heave  by  TRG.  This  method  uses  the  infinite- frequency 

potential  to  produce  a  simple  approximation  to  the  wave- amplitude 

coefficient  which  is  asymptotically  valid  both  at  low  and  high  frequencies. 

Numerical  calculations  in  a  number  of  cases  show  that  the  error  is  small 

at  all  frequencies.  (See  Figures  2,3). 

As  far  as  experimental  work  is  concerned,  very  little 

of  a  fundamental  nature  has  been  done  directly  on  the  two-dimensional 

r  13 1 

problem.  Stelson  has  measured  the  virtual  mass  of  a  heaving  sphere1  . 

An  earlier  set  of  investigations  of  both  virtual  mass  and  damping  of 
two-dimensional  forms  was  carried  out  by  Holstein but  there  was  a 
large  amount  of  scatter  in  the  data  so  that  comparisons  of  theory  and 
experiment  were  not  very  definitive.  What  appears  to  be  the  most 
useful  experimental  data  concerning  two- dimensional  virtual  mass  and 
damping  is  that  obtained  by  Porter*-  ■*.  The  total  hydrodynamic  force 
on  an  oscillating  semicircular  cylinder  was  obtained  and  compared  with 
theory,  showing  excellent  agreement.  While  no  separation  of  the  total 
force  into  virtual  mass  and  damping  components  was  possible  in  those 
tests  due  to  noise  and  other  instrumentation  imperfections,  it  appears 
that  this  result  adds  support  to  the  theoretical  results.  Nevertheless, 
a  precise  measurement  and  separation  of  the  virtual  mass  and  damping 
effects  for  such  a  fundamental  form  as  a  semicircle  would  serve  as  a 
basic  reference  for  comparison  with  theoretical  results,  and  it  is 
recommended  that  this  be  carried  out  . 


*  Our  attention  has  been  called  to  the  work  of  Yun-Sheng  Yu, 

Hydrodynamics  Laboratory,  M.I.T.,  who  compared  theory  and  experiment 
(wave  amplitude  and  added- mass)  for  a  heaving  semi-circle. 
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FIGURE  2 

RATIO  OF  WAVEHEIGHT  TO  HEAVE  AMPLITUDE 
FOR  ELLIPTIC  CYLINDERS.  H  =  HALF  -  BE  AM  /  DRAFT 
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2.2.2  Three-Dimensional  Methods.  Zero  Speed 

The  heaving  motion  of  a  semi- Immersed  sphere  has  been 
treated  by  HaveloclJ^,  Ursell^^,  Barakat^^,  TRG^J  and  porter^®^. 
Havelock  formulated  the  problem  as  an  Infinite  linear  system,  which  he 
truncated  and  solved  numerically.  Barakat  did  the  same  and  attempted 
to  perform  the  numerical  work  with  greater  accuracy  and  for  a  larger 
frequency  range.  Although  his  results  are  similar  to  Havelock's  there 
are  discrepancies,  and  there  Is  conclusive  evidence  In  Havelock's  favor. 
Ursell  found  the  high-frequency  behaviour  using  the  method  he  first 
applied  to  the  circular  cylinder,  and  found 


h(Havelock's  damping  parameter) 


added-mass  coefficient 


27 

4(Ka)4 

1 _ 3 

7  16ka 


Porter  has  very  recently  obtained  numerical  results  which  are  in 
excellent  agreement  with  Havelock.  Using  the  three-dimensional  form 
of  the  Hi-Fi  method  TRG  has  calculated  Havelock's  damping  parameter. 

The  results  are  shown  In  the  Figure  4. 

MacCamy  has  treated  the  heaving  circular  disc  by 
solving  an  Integral  equation  He  finds  the  wave  amplitude  and 

the  pressure  distribution  on  the  disc.  He  can  then  compute  the  damping 
coefficient  In  two  different  ways,  and  the  results  are  not  In  complete 
agreement.  TRG  has  computed  the  damping  coefficient  for  a  heaving 
circular  disc  by  the  Hi-Fi  method.  The  results  are  compared  in 
Figure  5. 

2.2.3  Three-Dimensional  Methods.  Forward  Speed 

The  problem  of  an  oscillating  surface  ship  at  zero 
speed  is  easily  linearized  by  assuming  either  a  thin  ship  or  small 
motions.  The  situation  with  forward  speed  is  more  complicated.  In 
the  method  of  Stoker  and  Peters,  £19],  who  perturb  the  nonlinear 
problem  for  small  beam,  various  interesting  quantities  are  zero  to 
first  order.  On  the  other  hand  a  perturbation  around  zero  speed  for 
a  ship  having  finite  dimensions  has  not  been  reported  in  the  literature. 
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FIGURE  4 

HAVELOCK'S  DAMPING  PARAMETER 
FOR  A  HEAVING  SPHERE 
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DAMPING  COEFFICIENT  FOR  A  HEAVING  CIRCULAR  DISK 
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For  a  submerged  body  linearization  of  the  free  surface  condition  has 
been  justified  on  the  basis  of  sufficient  distance  of  the  body  from 
the  free  surface,  together  with  small  motions.  However,  although  useful 
practical  results  have  been  obtained  in  this  manner  the  nature  of  the 
underlying  perturbation  theory,  if  any,  is  far  from  clear. 

Newman has  carried  out  a  perturbation  expansion 
in  which  the  amplitude  of  the  incident  waves,  the  amplitude  of  ship 
oscillation  and  the  ship1  s  beam  are  taken  as  independent  small  para¬ 
meters.  The  calculation  is  pushed  far  enough  to  yield  complicated  but 
nonvanishing  expressions  for  all  the  interesting  cruantities. 

A  technique  developed  by  Eggers^^  ^  gives  the  damping 
and  also  added  resistance  due  to  oscillatory  motions  of  sources  and/or 
dipoles  moving  with  forward  speed,  below  the  free  surface,  by  means  of 
an  energy  analysis.  The  problem  for  any  practical  application  is 
to  determine  the  necessary  distribution  strengths  and  to  carry  out  the 
required  integrations.  It  is  possible  that  Eggers'  work  will  lead  to 
a  useful  method  of  computing  the  da«pin|  of  submerged  bodies. 

More  recently  Vossers1  -1  has  published  a  s lender-body 
theory  (beam  and  draft  of  equal  small  order  compared  to  length)  in 
which  the  linearized  free  surface  condition  is  assumed  ab  initio. 

The  discussions  at  the  symposium,  on  this  subject, 
were  not  extensive  and  have  been  superceded  by  subsequent  work. 

As  far  as  experimental  results  are  concerned,  a  number 
of  investigations  using  oscillators  and  related  techniques  have  been 
carried  out  in  order  to  determine  the  virtual  mass  and/or  damping  of 
various  ship  hulls.  Haskind  and  Riman^^  carried  out  an  experimental 
study  which  resulted  in  values  of  virtual  mass  and  damping  of  a  simple 
ship  form  based  on  an  analytic  representation  as  a  member  of  a  simple 
family  of  ships.  The  results  obtained  therein  appeared  to  agree  quite 
closely  with  Haskind' s  earlier  theory,  but  there  is  insufficient 
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information  on  the  exact  data  as  well  as  the  possible  existence  of 
experimental  errors  which  are  not  considered  in  detail  in  that  study. 

A  later  set  of  experiments  on  a  Series  60  hull  was  carried  out  by 
Gerritsma^a*  ^  ^ ,  where  comparisons  were  made  with  the  results  of 
two-dimensional  theories  such  as  Havelock-Holstein  and  Grim,  which 
were  integrated  over  the  hull  in  a  strip  theory  application.  Better 
agreement  was  obtained  with  Grim's  results,  but  even  there  the  departure 
in  heave  damping  was  quite  large.  Experimental  measurements  on  a 
mathematical  ship  form  model  were  carried  out  with  an  oscillator  by 
Golcvato^®^  and  expressions  for  all  of  the  hydrodynamic  forces  and 
moments  due  to  heave  and  pitch  were  determined.  Newman^^  made 
comparisons  of  his  theoretical  results  for  damping  with  these  experi¬ 
mental  data  and  the  results  showed  large  differences,  with  the 
theoretical  results  much  higher  than  the  experiments.  While  some 
experimental  error  may  be  present,  the  order  of  the  differences 
indicates  that  a  fundamental  difference  is  evident.  Since  thin- ship 
theory  is  supposed  to  be  valid  only  for  thin  ships,  a  model  of  a 
thin  ship,  mathematically  defined  for  simplicity,  was  oscillated  by 
Gerritsma*^  ^  to  obtain  data  for  comparison  with  Newman's  theory. 

The  results  indicated  experimental  magnitudes  of  damping  larger  than 
theory,  and  led  to  speculation  about  the  relative  importance  of 
viscosity  affecting  measurements  of  damping  for  thin  forms. 

Most  of  the  work  on  seaworthiness  has  dealt  with  the 
vertical-plane  motions  of  heave  and  pitch  in  head  waves.  Similarly 
the  present  seminar  has  devoted  its  major  energies  toward  the  hydro- 
dynamic  theories  related  to  these  same  motions.  However,  it  appears 
that  lateral  motions  in  a  seaway,  which  are  related  to  course- keeping 
and  directional  stability1,  are  now  becoming  of  more  concern  and 
importance  as  the  vertical  plane  motions  are  understood  to  a  greater 
degree  than  initially.  Examining  the  available  literature  it  is  seen 

*  Recently,  Gerritsma,  Kerwin,  and  Newman  ^-1  used  the  same  theoreti- 
cal  and  experimental  techniques  to  determine  the  damping  of  the 
conventional  Series  60  forms,  and  in  this  study  the  agreement 
between  theory  and  experiments  was  very  much  improved. 
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that  two-dimensional  section  damping  and  virtual  mass  data  for  lateral 

T32  1 

motion  and  also  for  roll  are  available  in  the  work  of  GrimL  1  and  also 

T33 1 

in  some  Japanese  work1  .  The  importance  of  these  quantities  for  the 
case  of  zero  forward  speed  is  very  evident ,  but  serious  questions  arise 
as  to  their  relative  importance  at  speed.  This  is  due  to  the  occurrence 
of  large  hydrodynamic  fortes  and  moments  due  to  angle  of  attack  and 
angular  rotation,  which  arise  due  to  the  action  of  hydrodynamic  lift 
forces.  These  hydrodynamic  forces  are  much  larger  than  the  damping 
forces  due  to  oscillation  on  the  free  surface,  even  if  the  effect  of 
forward  speed  is  included  in  the  calculation  of  the  damping.  Thus 
emphasis  is  shifted  toward  the  dynamic  forces  that  determine  the  calm 
water  lateral  dynamic  stability  of  the  ship,  with  possible  small  modifi¬ 
cation  due  to  the  frequency  effects  (mainly  at  low  frequencies,  which 
are  the  most  important  for  lateral  motion) .  An  example  of  a  recent 
study  that  makes  use  of  some  of  these  ideas  is  that  by  Rydill^^,  but 
more  work  will  have  to  be  carried  out  in  order  to  put  this  approach 
on  a  more  rational  basis. 
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2.3  Exciting  Forces  on  Ships  in  Waves 

With  regard  to  the  problem  of  the  determination  of  the  exciting 
forces  acting  on  ships  in  waves,  the  major  effort  to  date  has  concerned 
itself  with  applying  the  results  obtained  for  submerged  bodies.  The 
application  of  slender  body  theory  to  the  evaluation  of  the  forces  acting 
on  submerged  bodies  under  waves  (neglecting  any  free  surface  effects 
and  self-wave  formation)  has  shown  good  success  in  achieving  a  simple 
representation  (see  the  work  of  Kaplan  and  Hu,  Recent  experimental 

studies ^23^  have  given  support  to  this  theory,  and  the  simple  interpreta¬ 
tion  in  terms  of  virtual  mass  and  fluid  acceleration  is  a  useful  tool 
in  further  extension  of  these  results.  Modification  of  the  virtual  mass 
term  to  account  for  frequency  dependence,  associated  with  free  surface 
influence,  leads  to  direct  application  to  the  surface  ship  problem,  as 
shown  by  Korvin-Kroukovsky  and  Jacobs ^  and  Kaplan^2^.  The  limited 
experimental  data  on  the  vertical  force  and  pitching  moment  on  a 
restrained  ship  model  in  regular  waves,  as  given  by  Korvin-Kroukovsky 
and  Jacobs,  shows  good  agreement  with  theory,  thereby  giving  support 
to  this  method.  It  is  intended  to  extend  this  method  to  other  forces, 
such  as  side  force  and  yawing  moment  in  oblique  waves,  with  some 
partial  account  of  diffraction  effects  in  allowing  the  ship  to  make  its 
own  waves.  Subsequent  to  the  seminar  Hu^33^  has  published  results  of 
this  type. 

As  far  as  three-dimensional  approaches  are  concerned,  the 
T251 

method  used  by  Cummins1  1  to  find  the  forces  on  a  submerged  body  under 
waves  is  presently  the  most  complete.  However,  a  knowledge  of  the 
singularity  distribution  representing  the  body  is  necessary,  and  there 
does  not  appear  to  be  any  simple  way  of  extending  his  method  to  the 
free  surface.  As  far  as  the  results  for  a  submerged  body  are  concerned 
however,  they  are  the  same  as  those  obtained  by  the  two-dimensional 
slender  body  technique  of  Kaplan  and  Hu1  J. 
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[97] 

Subsequent  to  the  Seminar  Newman1-  *  has  called  attention 
[281 

to  the  work  of  Haskind1-  who  showed  that  at  zero  speed  the  exciting 
forces  due  to  waves  can  be  determined  exactly,  including  the  diffraction 
effect,  without  solving  the  diffraction  problem.  In  fact  one  need  only 
know  the  wave  amplitude  coefficient  for  forced  motion  in  the  direction 
from  which  the  wave  is  coming  and  at  the  frequency  of  the  wave.  This  is 
a  fine  result  and  the  possibility  of  extension  to  forward  speed  is 
under  active  consideration  by  a  number  of  investigators. 
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2.4  Problems  and  RecommonMatjons 

2.4.1  Consistent  and  Inconsistent  Linearizations 

By  consistent  linearization  we  mean  the  systematic 
expansion  of  the  full  nonlinear  problem  In  terms  of  one  or  more  small 
parameters.  Examples  are  the  thin- ship  and  raft  theories  of  Stoker 
and  Peters,  and  Wehausen,  the  thin- ship  theory  of  Newman  and  a  theory 
for  slender  surface  ships  published  recently  by  Vossers.  Theories  In 
which  slender  submerged  bodies  are  treated  assuming  the  linearized  free 
surface  condition  ab  Initio  can  probably  be  made  consistent  If  we  regard 
them  as  resulting  from  an  expansion  In  two  small  parameters,  slenderness 
and  reciprocal  submergence,  although  the  nature  of  the  expansion  Is  not 
clear.  Inconsistent  theories  arise  whenever  we  assume  the  linearized 
free  surface  condition  ab  lnltlo  and  do  not  otherwise  restrict  the 
problem.  Much  of  the  work  on  submerged  finite  bodies  Is  of  this  type, 
e.g.  Havelock's  work  on  the  wave  resistance  of  a  submerged  circular 
cylinder,  MacCamy's  work  on  the  motions  of  shallow  cylinders  (rafts), 
strip  theory  based  on  section  coefficients  of  full  sections,  attempts 
at  computing  wave  resistance  assuming  the  linear  free  surface  condi¬ 
tions  while  satisfying  the  boundary  conditions  on  the  hull,  etc. 

It  was  the  feeling  of  the  participants  that  more 
problems  should  be  considered  by  Inconsistent  methods,  since  it  Is 
possible  thereby  to  get  concrete  results.  At  the  same  time  it  was 
suggested  that  the  time  has  coma  to  consider  the  simplest  cases  of 
a  body  (or  singularity)  using  the  exact  free  surface  condition,  or  a 
more  exact  approximation  to  it.  One  approach  suggested  was  to  consider 
a  line  dipole  (two-dimensional  problem)  numerically,  as  a  function  of 
the  dipole  strength  and  depth,  to  establish  the  range  of  strength-depth 
combinations  in  which  the  linearised  treatment  is  useful.  By  letting 
the  strltegth  increase  or  the  depth  decrease  one  would  hope  to  discover 
what  characteristically  nonlinear  features  appear.  It  was  also  proposed 
that  this  problem  be  treated  analytically. 
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2.4.2  Relations  Between  Cross-Coefficients  of  Inertia 
and  Damping 

In  the  treatment  of  a  body  which  is  at  the  same  time 
in  uniform  translation  and  undergoing  monochromatic  small  oscillations 
involving  six  degrees  of  freedom  two  six  by  six  matrices  arise.  One  of 
these  is  the  matrix  Bj.  of  damping* coefficients,  where  i  denotes  the 
degree  of  freedom  on  which  the  damping  force  acts  and  j  denotes  the 
degree  of  freedom  whose  velocity  gives  rise  to  the  damping  force.  The 
other  matrix  is  defined  similarly  but  involves  accelerations. 
Various  questions  arise  regarding  the  symmetry  properties  of  these 
matrices  and  the  vanishing  of  certain  elements  thereof.  The  dependence 
of  the  answers  on  the  type  of  linearization  assumed  and  on  the  symmetry 
properties  of  the  body  should  be  determined. 

After  the  conclusion  of  the  seminar,  Newman  and 
Tinman  have  looked  into  the  question  for  the  so-called  damping  matrix, 
assuming  that  the  body  has  transverse  and  longitudinal  symmetry  and  is 
sufficiently  thin,  slender  or  submerged  that  linear  theory  is  valid. 

The  results^2^  are  as  follows: 


Surge 

Sway 

Heave 

Roll 

Pitch 

Yaw 

>1 

j-2 

j-3 

j-4 

j-5 

j-6 

Surge  force 

i-1 

B11 

B12-0 

B13-B31 

b14-0 

B15*B51 

B16'° 

Sway  force 

i«2 

Heave  force 

B2i-° 

B22 

B23-° 

B24"B42 

B25*° 

B26"'B62 

1-3 

Roll  moment 

B3l"~B13 

B32"° 

B33 

B34-0 

B35-B53 

B36-0 

i-4 

B4l"° 

B42"B24 

*43-° 

B44 

B45’° 

B46"‘B64 

Pitch  moment 

i-5 

BsrBi5 

b52-° 

B53"'B35 

b54-0 

B55 

B56-° 

Yaw  moment 

i-6 

B6i-0 

B62"“B26 

B63*° 

B64""B46 

B65-° 

B66 

This  appellation  is  improper  when  there  is  coupling  between  two  or 
more  modes  of  oscillation. 
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Subsequently  Lurye^^,  has  investigated  the  symmetry  properties  of 
the  matrix  connecting  the  six  complex  forces  and  moments  with  the  six 
corresponding  complex  velocities. 

2.4.3  Relations  of  the  Kramers -Kronig  Type  Between  Damping 

and  Inertia  Coefficients 

At  the  seminar  it  was  suggested  by  Kotik  that  the  wave- 
amplitude  coefficient  and  the  added-mass,  as  functions  of  frequency,  are 
related  in  a  manner  similar  to  the  real  and  imaginary  parts  of  the 
complex  dielectric  constant  of  a  material,  or  the  resistive  and  reactive 
parts  of  the  impedance  of  a  linear  passive  network.  If  the  appropriate 
relation  were  found  it  would  be  very  useful.  For  instance,  it  would 
provide  an  internal  consistency  check  on  added-masses  and  wave- amplitude 
coefficients  computed  by  approximate  methods . 

Subsequent  to  the  seminar  the  details  were  worked 
cut  [30]  £0J.  thg  ca8e  of  zero  speed.  Some  numerical  consequences  are 
presented  in  Figures  6  and  7. 

2.4.4  Effect  of  Viscosity  in  Wave  Problems 

At  the  seminar  it  was  mentioned  that  certain  experi¬ 
ments  performed  by  Gerritsma  on  ship  forms  varying  in  beam  had  led  to 
peculiar  results  for  thin  models.  Some  aspects  of  the  experimental 
results  obtained  in  England  on  the  motions  of  thin  planks  are  also 
not  predicted  by  theory.  It  was  suggested  that  viscosity  might  be 
important  for  thin  models,  and  although  no  conclusion  was  reached  the 
following  theoretical  problems  were  proposed;  to  find  the  waves 
generated  by  a  heaving  vertical  plate  of  zero  thickness,  and  by  a 
rolling  circular  cylinder.  These  problems  are  of  special  interest 
because  in  the  absence  of  viscosity  no  waves  would  be  produced.  After 
the  seminar  Karp  found  that  an  approximate  expression  for  the  wave 
amplitude  coefficient  for  a  heaving  vertical  plate  in  a  viscous  fluid, 
valid  when  k /  «  1,  is 
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^  TRG,  OBTAINED  BY  SUBSTITUTING  THE  HI-FI 
DAMPING  COEFFICIENT  INTO  THE  KRAMERS- 
KRONIG  RELATION. 


FIGURE  6 

ADDED- MASS  PARAMETER  FOR  A 
HEAVING  CIRCULAR  CYLINDER 


6  0 


FIGURE  7 

ADDED-MASS  PARAMETER  FOR  A  HEAVING  SPHERE 
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2k/, 

where 

k  -  cr2/g 

and 

/  ■  -Jv/a 

is  a  measure  of  the  boundary  layer  thickness  in  an  infinite  fluid. 

2.4.5  Time-Dom<Hn  Formulation  of  the  Equations  of  Motion 
At  the  seminar  the  question  of  the  form  of  the 

equations  of  motion  for  nonmonochromatic  excitation  was  brought  up,  and 
the  usual  remarks  regarding  the  frequency- dependence  of  the  damping  and 
added-mass  coefficients  were  made.  It  was  suggested  that  the  appropriate 
equations  would  be  integro-differential  equations  with  convolution 
kernels.  Subsequently  Cummins has  supplied  the  appropriate  analysis, 
and  the  questions  have  been  largely  resolved,  in  the  linear 

range  the  time  and  frequency- domain  formulations  are  completely 
equivalent . 

2.4.6  Rec  ommendat ions 

It  appears  that  a  number  of  the  problems  put  forth 
at  the  seminar  have  subsequently  been  solved,  and  the  seminar  may  be 
credited  with  a  directly  stimulating  role  in  this  regard.  Most  of  the 
participants  felt  that  they  benefited  in  various  ways  not  concretely 
reflected  in  subsequent  publications  but  valuable  none  the  less.  It 
appears  that  a  small  meeting  of  specialists,  like  this  seminar,  can 
fruitfully  be  held  at  Intervals  of  about  three  years. 

The  value  of  such  a  seminar  to  the  field  as  a  whole 
could  be  improved  by  introducing  a  suitable  mechanism  for  recording 
the  principal  points  brought  out  in  discussion.  The  method  used  in 
this  seminar  was  to  tape  record  the  discussion.  Analysis  of  the  records 
was  unfortunately  delayed,  and  the  method  is  in  any  case  not  satisfactory. 
It  is  recommended  that  in  the  future  one  of  the  participants  be  assigned 
to  write  up  each  half-day  of  discussion,  immediately  after  it  occurs. 

The  writeups  should  be  circularized  and  amended  during  the  seminar  and 
before  it  concludes,  and  the  final  version  should  be  made  available  for 
distribution  shortly  thereafter. 
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